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Abstract. Diagonal and off-diagonal flavor and conserved charge susceptibilities 
reveal the prevalent degrees of freedom of heated strongly interacting matter. 
Results obtained from lattice simulations are compared with various model 
estimates in an effort to weed down various possible pictures of a quark 
gluon plasma. We argue that the vanishing of the off-diagonal quark flavor 
susceptibilities and its derivatives with respect to chemical potential, at 
temperatures above 1.5 Tc, can only be understood in a picture of a gas or 
liquid composed of quasi-particles which carry the same quantum numbers as 
quarks and antiquarks. A potential new observable, blind to neutral and non- 
strange particles, is introduced and related via isospin symmetry to the ratio of 
susceptibilities of baryonic strangeness to strangeness generated in the excited 
matter created at RHIC. 



PACS numbers: 12.38.Mh, 11.10. Wx, 25.75.Dw 

The goal of of the heavy-ion program at the Relativistic Heavy Ion Collidcr 
(RHIC) is the creation and study of heated strongly interacting matter at nearly 
vanishing baryon density [1, 2]. The experimental results of the four RHIC detector 
collaborations, have set a lower bound of about 5 GeV/fm 3 on the energy density 
at a time r = 1 fm/c in central Au+Au collisions [3]. Such an energy density 
should place the produced matter well into region of the QCD phase diagram that 
cannot be described as a hadronic resonance gas. Matter in this domain of the QCD 
phase diagram had been expected to be a colored plasma composed of quasiparticle 
excitations with the quantum numbers of quarks and gluons [1, 4]. The produced 
matter, however, clearly exhibits collective behavior as evidenced by its radial and 
elliptic flow [5] . This result suggests that there must exist a strong interaction between 
the constituents of the medium. 

Lattice simulations of QCD at finite temperature and vanishing chemical potential 
have demonstrated the existence of a phase transition at T = T c ~ 170McV signalled 
by a steep rise in the energy density and the pressure as a function of the temperature. 
The slow rise of both quantities prior to the sudden transition has been successfully 
modelled by a picture of a hadronic resonance gas [6]. A similar program for the 
description of the excited phase as a weakly interacting plasma of quasiparticles [7] 
has not met with success in the region from T c < T < 3T C . This indicates that 
matter in this region may not be a weakly coupled plasma where quarks and gluons 
arc deconfined over large distances. It is clear that a microscopic understanding of 
the emergent degrees of freedom in this regime is essential in order to understand 
the matter created in nuclear collisions at RHIC. The object of these proceedings 
is to outline a set of flavor off-diagonal susceptibilities which may be calculated 
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on the lattice, in phcnomcnological models and in certain situations may even be 
measured in experiments. Such off-diagonal susceptibilities are discerning probes of 
the thermodynamic degrees of freedom in strongly interacting QCD matter, both 
above and below T c . 

In the lattice formulation of QCD, the fundamental degrees of freedom are local 
quark and gluon fields. Under conditions where deconfinement has been achieved, the 
elementary set of conserved charges is given by the net contents of each quark flavor 
u,d,s. An alternate basis is provided by the hadronically defined conserved charges 
of baryon number (B), electric charge (Q) and strangeness (S). The (co-)variances 
between any two conserved charges (x, y) are extensive quantities: a\ y = VTxxy, 
where Xxy is the intensive diagonal or off-diagonal susceptibility. These susceptibilities 
can be measured on the lattice. In heavy-ion experiments, the variances and 
covariances are measured by means of an event- by- event analysis of the corresponding 
conserved quantities, i.e., 

where E represents the set of events, Ne is the number of events considered, and 
Xi,Yi are the net values of the conserved charge in a given event i. The volume 
independent ratios of variances, measured event-by-event in heavy-ion collisions, may 
then be directly compared with the lattice estimate for the ratio of susceptibilities. 

In the case that the degrees of freedom or active flavors / are eigenstates of the 
quantum numbers x, y, the above expression may be decomposed using Xi — njx/ 

where n{ are the number of states of flavor / in event i. For independent flavors, where 
(n^n 9 ) = J2ieE n i n i = (nf)(n 9 ), one may easily demonstrate that (sec Rcf. [8] for 
details) 



The last equality in the above equation holds solely in the case that Poisson statistics 
is applicable to the independent flavors (i.e., the variance a 1 is equal to the mean (n)) 
Equation (2) clearly demonstrates the sensitivity of the off-diagonal susceptibilities to 
the presence of flavors carrying more than one quantum number, and may be used to 
estimate the susceptibilities in various microscopic models. 

A robust set of ratios of susceptibilities, which has been measured on the lattice 
and demonstrates direct sensitivity to the flavor carrying degrees of freedom are 
the two coefficients (see Refs. [9, 10]) Cbs = —3°~Bs/ a s an d Cqs = 3°~Qs/ a s- 
The susceptibilities may be computed phenomenologically, assuming the condition 
of Eqs. (2). Plotted in the left panel of Fig. 1. are the coefficients Cbs an( i Cqs as 
obtained from a hadronic resonance gas spectrum, truncated at the mass of the Q~ 
along with the results obtained from the lattice. One notes, that the hadron resonance 
gas provides a good description of the behavior of the ratio of susceptibilities up to 
the point of the phase transition. Here the behavior of the truncated spectrum fails 
to reproduce the sharp rise in Cbs and the corresponding sharp drop in Cqs- 

The results from the lattice simplify in the high temperature phase, where one 
notes the following simple relations [11, 12]; 




(2) 



Xus = Xds < Xud < Xs < Xd = Xu, 



(3) 
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Figure 1. (Left Panel) A comparison of CbSi Cqs calculated in a truncated 
hadron resonance gas at /is = us = fiQ = OMeV compared to lattice calculations 
from Ref. [10]. The two hazed bands for Cbs an d Cqs for the hadron gas plots 
reflect the uncertainty in the actual value of the phase transition temperature 
T c = 170 ± lOMeV. 

(Right Panel) A plot of the mesonic and baryonic contributions to Xud{T,p = 0) 
and and a plot of the derivatives Td 2 Xud/9(fJ,/3) 2 \ l _ l= o (which is equal to the 
baryonic contribution pg) and Td 4 Xud/d(l J -P) 4! \^=0 (which is negative beyond 
T c = 170MeV and thus is inconsistent with a bound state interpretation). 



where the last and the first equality applies in the case of isospin symmetry. For 
simulations at vanishing chemical potentials (B) = (S) = (Q) = 0, as a result, we 
obtain, 

n (BS) ((u + d + s)s) Xus + Xds + Xs , , A s 
Cbs - -3 j-wr 1. (4) 

Alternatively, one may state that Xus — Xds — * 0. Models of the deconfined phase 
must obey such constraints. The simplest model of deconfined matter is that of 
non-interacting quark, anti-quark and gluon quasi-particles. In such a situation, off- 
diagonal susceptibilities are identically zero as the degrees of freedom may carry only 
a single flavor. 

Recently, an alternate picture of the deconfined phase has been proposed, where 
a tower of bound states of quarks and gluons arc present besides the quasiparticles 
themselves [13]. The large cross-sections from resonance scattering imply short mean 
free paths and lend consistency to the macroscopic hydrodynamic picture henceforth 
used to describe the dynamical evolution of the matter produced at RHIC. However, 
in a plasma containing bound states of quarks, in the form of colored mesons, diquarks 
and quark-gluon bound states, the correlation between flavors is no longer negligible. 
In a gas consisting solely of quark-antiquark bound states, built from three flavors 
of quarks, the contributing states are ud, du, us, su, ds, sd (The flavor singlets qq and 
gluon bound states will be ignored as they carry no conserved flavor). In such a 
system, the ratio Cbs must vanish, because all states have vanishing baryon number. 
Thus meson like states produce a considerable negative contribution to the off-diagonal 
susceptibilities Xus,Xud,Xds and tend to depress the value of Cbs from unity. Quark 
gluon bound states contribute similarly as quark-antiquark quasiparticles. Inclusion 
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of all such states at a temperature of T = 1.5T C yields a value Cbs — 0.62 (or 
equivalcntly, negative off-diagonal susceptibilities) quite different from the values 
found on the lattice [9] . 

These results have motivated the inclusion of a variety of baryonic states into 
the model outlined above [14]. Baryonic states provide additive contributions to the 
off-diagonal susceptibilities and as a result also to Cbs- The inclusion of a sufficient 
number of baryonic states may lead to a balance between the mesonic and baryonic 
contributions and one may engineer a Xus = Xds = and as a result a Cbs = 1 
in agreement with the lattice. To discern between these two possibilities, one turns 
to higher derivatives of the the off-diagonal susceptibilities in terms of the baryon- 
chemical potential ps- 

In the interest of simplicity, the flavor group will be restricted to SU(2)/. Lattice 
results for susceptibilities and their derivatives with respect to baryon chemical 
potential using dynamical quarks exist in this case [12]. The model will consist of quark 
and anti-quark quasi-particles u, d and u, d, meson like bound states uu, dd, ud, du, 
diquark states uu, dd, ud and their antiparticles as well as baryons uuu, uud, udd, ddd 
and the corresponding anti-baryons. In this simplified situation, one may derive the 
following simple relations between the off-diagonal susceptibility and the meson and 
baryon densities [8] (p% = n®/V, where n° denotes the population of a given flavor at 
1-1 = and not the difference between the flavor anti-flavor populations) 

TXud(T^)\ =0 ~ - 2p° ud (T) + {2p° ud (T) + i P ° uud (T) + ip° udd (T)}, (5) 



T 



d 2 Xud 



d{ppy 



d 4 Xua 



(6) 



= 2p° ud (T) + Apl ud {T) + 4p° dd (T) = T 

ju=0 

In this way, one may divide the contributions to the off-diagonal susceptibility 
and its derivatives in terms of mesonic [pm = 2p d (T)) and baryonic (ps = 
^Pud( T ) + ^Puud( T ) + 4Pudd( T )) contributions. Using the measured susceptibility 
and its derivatives, these densities are estimated as a function of the temperature, 
and plotted in the right panel of Fig. 1. These densities satisfy Eq. (5) and the first 
equality of Eq. (6). The second condition imposed by Eq. (6) has to be satisfied by the 
fourth derivative of the susceptibility in such a picture of bound states. The fourth 
derivative of the susceptibility, as measured on the lattice, has been plotted as the 
square symbols. Despite large error bars, one notes that while the baryon density or 
the second derivative of the susceptibility is consistent with Eq. (6) below the phase 
transition temperature, it becomes inconsistent with such a condition above the phase 
transition temperature. Above T c , Td 4 Xud/d(p/3) 4 is actually negative; hence, no 
bound state picture is compatible with such results. It has already been pointed out 
in Rcf. [12] that the signs of the various derivatives of the susceptibility arc consistent 
with the picture of a weakly interacting quasi-particlc gas. 

In the experimental measurement of baryon-strangeness correlations, any detector 
has to accurately assess the baryon number and strangeness in a given rapidity bin 
in each event. As most detectors are blind to stable uncharged particles they cannot 
measure the neutron and antineutron populations. As a result, a measurement of o~bs 
may become rather difficult. As a recourse, a new quantum number is constructed, 
M = B + 2I3, and fluctuations of M with respect to S are studied. One notes that the 
assumption of isospin symmetry reduces the covariance ctms to simply o~bs (see Ref. [8] 
for details). As a result, in all theoretical models with isospin symmetry Cms = Cbs- 
In the experimental determination, M has the advantage that it is vanishing for 
all particles that do not carry charge or strangeness, thus M = for neutrons, 
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Figure 2. A comparison of the two related ratios of variances Cbs, Cms as a 
function of the acceptance in rapidity from —\y m ax\ to |j/max|- 



antincutrons, neutral pions, etc.. There exists a greater possibility of contamination 
of such a signal from charged pions (tt^ do not contribute to the measurement of 
Cbs)- The amount of contamination is assessed by comparing measurements of Cbs 
and Cms m an identical set of HIJING events [15]. 

In Fig. 2 the correlations Cbs an d Cms are estimated in central Au — Au events 
at ^/s = 200 AGcV. The acceptance in rapidity ranges from — |y ma a;| to |y m aa;|, hence 
a larger y max indicates a larger acceptance. The results are presented as a function of 
Umax- One notes that over the range of y m ax the two correlations Cbs and Cms are 
rather similar. This bodes well for the measurement of Cbs m RHIC experiments via 
a measurement of the quantity Cms> over a range of rapidity intervals. 

In these proceedings, we have demonstrated that the off-diagonal susceptibilities 
such as Xud, Xus,Xds and the ratio Cbs are descerning probes of the flavor structure of 
heated strongly interacting matter. They may be calculated on the lattice, in models 
and even measured in RHIC experiments allowing for a systematic exploration of 
degrees of freedom of deconhned matter. 
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